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Q ■ Abstract 

Q^' Using results on soft-coUinear factorization for inclusive i?-meson decay distributions, 

a systematic study of the partial B Xgj decay rate with a cut > Eq on photon 
D ! energy is performed. For values of £"0 below about 1.9 GeV, the rate can be calculated 

' without reference to shape functions using a multi-scale operator product expansion 

(MSOPE). The transition from the shape-function region to the MSOPE region is stud- 
^ , led analytically. The resulting prediction for the B Xg'j branching ratio depends 

I on three large scales: rrib, \/m;,A, and A = nih — 2Eq. Logarithms associated with 

these scales are resummed at next-to- next-to-leading logarithmic order. While power 
corrections in Aqcd/^ turn out to be small, the sensitivity to the scale A ~ 1.1 GeV 
(for Eq ~ 1.8 GeV) introduces significant perturbative uncertainties, which so far have 
been ignored. The new theoretical prediction for the B Xg^ branching ratio with 
E^ > 1.8 GeV is Br(5 ^ ^^7) = (3.38tg-j^;[!;^g) x 10"'', where the first error is an 
estimate of perturbative uncertainties and the second one reflects uncertainties in input 
parameters. With this cut (89 ±1)% of all events are contained. The implications of 
larger theory uncertainties for New Physics searches are briefly explored with the example 
of the type-II two-Higgs-doublet model, for which the lower bound on the charged-Higgs 
mass is reduced compared with previous estimates to approximately 200 GeV at 
confidence level. 



1 Introduction 



The inclusive, weak radiative decay B — Xg'y is the prototype of all flavor- changing neutral 
current processes. In the Standard Model, this process is mediated by loop diagrams containing 
W bosons and top (or lighter) quarks. In extensions of the Standard Model, other heavy 
particles propagating in loops can give sizable contributions, which in many cases can compete 
with those of the Standard Model. As a result, measurements of the B Xs'~f rate and 
CP asymmetry provide sensitive probes for New Physics at the TeV scale. In many cases, 
the fact that these measurements agree with Standard Model predictions imposes non-trivial 
constraints on the allowed parameter space. 

Given the prominent role of i? — > Xg^ decay in searching for physics beyond the Standard 
Model, it is of great importance to have a precise prediction for its inclusive rate and CP 
asymmetry in the Standard Model. This has been achieved thanks to the combined effort of 
many theorists over a period of several years [1]. The total inclusive rate is known at next-to- 
leading order in renormalization-group (RG) improved perturbation theory with a theoretical 
precision of about 10%. Currently, a major effort is underway to improve this accuracy by 
calculating the dominant parts of the next-to- next-to-leading corrections [2, 3]. 

While the total inclusive B Xgj decay rate can be calculated using a conventional 
operator-product expansion (OPE) based on an expansion in logarithms and inverse powers 
of the 6-quark mass [4], the situation is more complicated when a cut on the photon energy 
is applied. In practice, experiments can only measure the high-energy part of the photon 
spectrum, > Eq, where typically Eq = 2 GeV (measured in the B-meson rest frame) or 
slightly lower [5, 6]. Even if such a cut was not required for experimental reasons, it would 
be needed to reduce the photon background from B — > Xgip^'^ decays followed by a radiative 
decay of the ■0^'^ [7]. With Ej restricted to be close to the kinematic endpoint at Mb/2 
(neglecting the kaon mass), the hadronic final state Xs is constrained to have large energy 
Ex ~ Mb but only moderate invariant mass Mx ~ (MbAqcd)^^^- In this kinematic region, 
important hadronic effects need to be taken into account. An infinite set of leading-twist terms 
in the OPE need to be resummed into a non-perturbative shape function, which describes the 
momentum distribution of the 6-quark inside the B meson [8, 9, 10]. In addition, Sudakov 
double logarithms arise near the endpoint of the photon spectrum, which need to be resummed 
to all orders in perturbation theory [11, 12, 13]. While these issues are now well understood 
theoretically [14, 15], the presence of the shape function leads to an unavoidable element of 
hadronic uncertainty and modeling, which is undesirable when the goal is to probe for physics 
beyond the Standard Model. 

Conventional wisdom says that, while shape- function effects are important near the end- 
point of the photon spectrum, these effects can be ignored as soon as the cutoff Eq is lowered 
below about 1.9 GeV. This assumption is based on phenomenological studies of shape-function 
effects using various model functions, which have the unrealistic feature that the distribution 
function vanishes exponentially for large light-cone momenta [16, 17]. In other words, it has 
been implicitly assumed that there is an instantaneous transition from the "shape-function 
region" of large non-perturbative corrections to the "OPE region", in which hadronic cor- 
rections to the rate are suppressed by at least two powers of Aqco/'^ife- As a result, the 
preferred strategy has been to encourage experimenters to lower the photon-energy cut to a 
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value E'o < 1.9 GeV, and then to employ the conventional OPE for the calculation of the rate, 
ignoring shape-function effects. 

In this paper we show that this strategy is based on a misconception. Our work is motivated 
by two considerations. First, is has recently been shown that the asymptotic behavior of B- 
meson distribution functions such as the shape function is not exponential, but rather governed 
by radiative tails exhibiting a slow, power-like fall-off [14, 18]. One should therefore not exclude 
the possibility of a significant radiation tail in the case of the B Xg'-f photon spectrum, 
meaning that more events than predicted by existing models could be located at low photon 
energy. Fits to experimental data in the low-energy part of the spectrum, which are based 
on such models, should thus be taken with caution. Secondly, it has been our desire for a 
long time to find an analytic way to study the transition from the shape-function region to 
the OPE region. If it were true that shape-function effects become irrelevant once the cutoff 
Eq is lowered below 1.9 GeV, one should be able to see this analytically using some form of 
a short- distance expansion. We show that this expansion indeed exists, and that it involves 
three different short-distance scales. In addition to the hard scale rri},, an intermediate "hard- 
coUinear" scale \/rr%K corresponding to the typical invariant mass of the hadronic final state 
Xsi and a low scale A = — 2Eq related to the width of the energy window over which the 
measurement is performed, become of crucial importance. The physics associated with these 
scales can be disentangled using recent results on soft-collinear factorization theorems derived 
in the framework of effective field theory [14, 15]. A systematic treatment consists of matching 
QCD onto soft-collinear effective theory (SCET) [19] in a first step, in which hard quantum 
fiuctuations are integrated out. In a second step, hard-collinear modes are integrated out by 
matching SCET onto heavy-quark effective theory (HQET) [20]. Ultimately, the precision of 
the theoretical calculations is determined by the value of the lowest short-distance scale A, 
which in practice is of order 1 GeV or only slightly larger. The theoretical accuracy that can 
be reached is therefore not as good as in the case of a conventional heavy-quark expansion 
applied to the B system, but more likely it is similar to (if not worse than) the accuracy 
reached in the description of the inclusive hadronic r decay rate i?^ [21]. However, while the 
ratio Rt is known to order a^, the B — > branching ratio is currently only known through 
order ctg. 

While we are aware that this conclusion may come as a surprise to many practitioners in 
the field of flavor physics, we believe that it is an unavoidable consequence of the analysis 
presented in this paper. Not surprisingly, then, we find that the error estimates for the partial 
B — > Xg'^ branching ratio in the literature are too optimistic. Since there are unknown ^^(A) 
corrections at the low scale A ~ 1 GeV, we estimate the present perturbative uncertainty in the 
B — > Xg^ branching ratio with Eq in the range between 1.6 and 1.8 GeV to be of order 10%. In 
addition, there are uncertainties due to other sources, such as the h- and c-quark masses. The 
combined theoretical uncertainty is of order 15%. While this is a rather pessimistic conclusion, 
we stress that the uncertainty is limited by unknown, higher-order perturbative terms, not 
by non-perturbative effects, which we find to be under good control. (This is similar to the 
case of Therefore, there is room for a reduction of the error by means of well-controlled 
perturbative calculations. 

In Section 2, we discuss the QCD factorization formula for the partial B X^'-y decay rate 
with a cut E^ > Eq on photon energy, valid at leading power in the heavy-quark expansion. 
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Contributions associated with the hard, hard-coUinear, and soft scales are factorized into 
a hard function a jet function J, and a shape function S. Single and double (Sudakov) 
logarithms are systematically resummed to all orders in perturbation theory. The RG evolution 
of the shape function is studied in Section 3, where we present the exact solution to its 
evolution equation in momentum space. Our main results are derived in Section 4, where we 
show how the convolution integral over the shape function in the factorization formula can 
be calculated using a local OPE, provided that the scale A = mi — 2Eq is numerically large 
compared with Aqcd- Section 5 discusses how to eliminate the HQET parameters rrif, and 
Ai defined in the pole scheme, which enter the theoretical expressions, in terms of physical 
parameters defined in the so-called "shape-function scheme" [14] . The calculation of the decay 
rate is completed in Section 6, where we add contributions that are power-suppressed in the 
heavy-quark expansion. For these small corrections, the scale separation we achieve is only 
approximate and misses some yet unknown terms of order ln^(A/m;,). In Section 7, we 
show that by considering ratios of decay rates one may separate the short-distance physics 
contained in the hard function from the physics associated with the intermediate and low 
scales. For instance, at leading power in A/m^ the ratio of the B — > Xg^ branching ratio 
in a New Physics model relative to that in the Standard Model can be calculated without 
any sensitivity to scales less than the hard scale fi^ ~ "^6, and the same is true for the 
direct CP asymmetry. In other cases, at leading power some ratios are insensitive to the 
Wilson coefficients in the effective weak Hamiltonian and thus to New Physics. Examples are 
the average photon energy (Ej), and the ratio of the B — > Xg^ decay rate with Ej > Eq 
normalized to the total rate. The latter ratio is particularly interesting, since it can be used 
to make contact between a simple, fully inclusive rate calculation and our more sophisticated 
analysis of multi-scale effects. Our numerical results are presented in Section 8, followed by a 
summary and conclusions. 

2 QCD factorization theorem 

Recent results on the factorization of hard, hard-collinear, and soft contributions to inclusive 
S-meson decay distributions [14, 15] allow us to obtain a QCD factorization formula for the 
integrated B Xs^ decay rate with a cut E^^ > Eq on photon energy. In the region of large 
Eq, the leading contribution to the rate can be written in the form 

X / dP+ {Mb — P+Y / duj nib J{''^h{P+ — '^), A*i) S{uj, /li) + power corrections, 
Jo Jo 

where Ag = Mb — 2Eq is twice the width of the window in photon energy over which the 
measurement of the decay rate is performed. In the prefactor, a is the fine-structure constant 
normalized a,t — [22]. The variable P+ = Ex — \Px\ is the "plus component" of the 
4- momentum of the hadronic final state Xg, which is related to the photon energy by P+ = 
Mb — 2E^. The factor {Mb — P+)^ under the integral thus equals 8E^, where two powers of 
E^ come from the squared matrix element of the effective weak Hamiltonian, and one factor 
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comes from phase space. The hadronic invariant mass of the final state is M|- = MbP+- The 
endpoint region of the photon spectrum is defined by the requirement that < ^ Mb, 
in which case is called a hard-collinear momentum [23]. Power corrections to the expression 
above will be analyzed later; however, the leading non-perturbative corrections to the total 
decay rate have already been factored out in (1) and included in the parameter [4] (see Table 1 
in Section 8 for a list of input parameters) 

(2) 

The factorization formula (1) was first presented in [12]. What is new is that we now have a 
systematic effective field-theory technology to compute the functions and J order by order 
in perturbation theory, and to control their scale dependence in momentum space (not moment 
space). Also, it is in principle possible to include power-suppressed terms in the heavy-quark 
expansion. In the factorization formula, fi^ ~ mi, is a hard scale, while /Xj ~ a/'"^6Aqcd is an 
intermediate hard-collinear scale of order the invariant mass of the hadronic final state. The 
precise values of these matching scales are irrelevant, since the rate is formally independent of 
Hh and The hard corrections captured by the function H^{f^h) result from the matching 
of the effective weak Hamiltonian of the Standard Model (or any of its extensions) onto a 
leading-order current operator of SCET. It is defined by the relation 
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^ ^tbVtl ^ E^m,{i,H) HM e;{q) [^WhcY±{'^ - 75) K] ilih) + ■ ■ ■ , (3) 

where e{q) is the transverse photon polarization vector, and the dots represent power-sup- 
pressed contributions from higher-dimensional SCET operators. The result is proportional to 
the photon energy, — v • q, defined in the S-meson rest frame. (Here v is the 4-velocity 
of the B meson.) The fields and ^ represent the soft heavy quark and the hard-collinear 
strange quark, respectively, and Whc is a Wilson line. At tree level, only the dipole operator 
Qr^ and the four-quark penguin operators Q5 and Qe in the effective weak Hamiltonian give a 
non-zero contribution to H^, which is equal to the "effective" coefficient Cy^ = — ^C^ — Cq. 
(We use the conventions of [24] for the operators and Wilson coefficients in the effective weak 
Hamiltonian.) At next-to-leading order, the result reads (with Cp — 4/3) 



1 H -2 In h 71n 6 - — ] + e^w 

47r \ Hh IJ'h 12/ 
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The coefficient Cj^{iJ,h) of the electromagnetic dipole operator is required with next-to- leading 

order accuracy [25], while the remaining coefficients can be calculated at leading logarithmic 
order. Explicit expressions for these coefficients can be found, e.g., in [17, 25]. The terms in the 
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third row arise from charm-quark and up-quark penguin contractions of the current-current 
operators Qi These contributions depend on the small ratio 
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X\p-if)) 1 + X\l - p - if]) + 0{X^ 



(5) 



The variable z = {rric/mij)'^ denotes the ratio of quark masses relevant to the charm loop, and 
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-9 - 147r^ + 1127ri + (182 - 487ri) In z - 126 In^ z 



+ 



(6) 



are the first few terms in the expansion of the penguin function [26], whose exact expression 
(in the form of parameter integrals) can be found in [27]. The imaginary parts in (4) and (6) 
are strong-interaction phases, which in conjunction with CP-violating weak phases contained 
in the parameter £ckm or in potential New Physics contributions to the Wilson coefficients 
can induce a non-zero CP asymmetry in S — > Xg^ decays [28, 29]. The term 



-1. 



(7) 



accounts for electroweak matching corrections at the weak scale [30] and logarithmically en- 
hanced electromagnetic effects affecting the evolution of the Wilson coefficients [17, 31]. Fi- 
nally, the term epeng ~ 0.2% includes the effects of penguin contractions of operators other 
than [27], which are numerically neghgible but are included here for completeness. In the 
factorization formula (1), the hard function is multiplied by the running 6-quark mass 



rribiiih) = mb{mb) 



1 + 



27r 



ln^ + 
Ph 



(8) 



defined in the MS scheme, which is part of the electromagnetic dipolc operator Qj-^. On the 
other hand, the scheme to be used for the quark masses entering the ratio z in the penguin 
function g{z) is not specified at next-to-leading order [7]. Since the matching is performed at a 
hard scale the charm-quark mass should be a running mass mc{ph)-i while enters either 
as the mass in the 6-quark propagator or via the values of external momenta. For simplicity, 
we take z = \mc{ph)) /^b{ph)]'^ as a ratio of running quark masses evaluated at the same scale, 
which has the advantage that this quantity is RG invariant. 
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The jet function J{mb{P+ — '^), A*i) in (1) describes the physics of the final-state hadronic 
jet. At next-to-leading order in perturbation theory, it is given by the expression [14, 15] 



rub J{mb{P+ - u}),i^i) = 6{P+ - a)) 



1 + 
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An 



P+-Cj 



[Mi Mi.] 



(9) 



The star distributions are generalized plus distributions defined as [32] 

1 
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(10) 



where F[x) is a smooth test function. The pcrturbative expansion of the jet function can be 
trusted as long as /if ~ m^A with A ~ pmax _ ^/^^ ^ _ -c Mb- By quark-hadron 
duality, only the maximum values of kinematic variables such as P+, which are integrated over 
phase space, matter for the calculation of inclusive decay rates [4] . Note that the ''natural" 
choices oc rUb and ^1 = nib P'i with /ij independent of remove all reference to the 6-quark 
mass (other than in the arguments of running coupling constants) from the factorization 
formula (1). 

The shape function S'(a), Hi) parameterizes our ignorance about the soft physics associated 
with bound-state effects inside the B meson [8, 9]. Its naive interpretation is that of a parton 
distribution function, governing the distribution of the light-cone component fc+ of the residual 
momentum k = Pb — rUbV of the b quark inside the heavy meson. Once radiative corrections 
are included, however, a probabilistic interpretation of the shape function breaks down [14]. 
For convenience, the shape function is renormalized in (1) at the intermediate hard-collinear 
scale rather than at a hadronic scale /^had- This removes any uncertainties related to the 
evolution from /li to /ihad- Since the shape function is universal, all that matters is that it is 
renormalized at the same scale when comparing different processes. 

The last ingredient in the factorization formula (1) is the RG evolution function f/i(/x/j, /Xj), 
which describes the evolution of the hard function \Hy\'^ from the high matching scale fi^ down 
to the intermediate scale /Xj, at which the jet and shape functions are renormalized. The exact 
expression for this quantity follows from 



rrib 



\nUi{i^h: l^i) = 2S{iih:l^i) - 2ar(/i/i,/ii)ln 2ay(/i/j,/i, 



(11) 



where the various functions on the right-hand side are the solutions to the partial differential 
equations (the Sudakov exponent S should not be confused with the shape function S) 

S{iy,ii) = -rcusp(tts(/^)) In-, 
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a^>{u,jj) = -7'(as(/^)) , 



(12) 
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with initial conditions S{i', u) = ar(i^, i') = ay(i/, i/) = at /i — u. Here Fcusp is the universal 
cusp anomalous dimension for Wilson loops with light-like segments [33], which has recently 
been calculated to three-loop order [34], and 7' enters the anomalous dimension of the leading- 
order SCET current operators containing a heavy quark and a hard-coUinear quark with large 
energy E, which takes the form [19, 23] 



As explained in Appendix A.l, a conjecture for the two-loop expression for 7' can be de- 
duced using results from the hterature on deep- inelastic scattering [35, 36, 37]. The evolu- 
tion equations (12) are solved in the standard way by writing d/d\iafi = (5{as) d/dag, where 
P{as) — das /din /I is the QCD /? function. This yields the exact solutions 



and similarly for the function ay. The perturbativc expansions of the anomalous dimensions 
and the resulting expressions for the evolution functions valid through order ag are collected 
in Appendix A.l. 

As written in (1), the decay rate is sensitive to non-perturbative hadronic physics via 
its dependence on the shape function. This sensitivity is unavoidable as long as the scale 
A = mfe — 2_Eo is a hadronic scale, corresponding to the endpoint region of the photon spectrum 
above, say, 2 GeV. The properties of the B Xg'y decay rate and photon spectrum in this 
region will be discussed in detail elsewhere. Here we are interested in a situation where Eq 
is lowered out of the shape-function region, such that A can be considered large compared 
with Aqcd- For orientation, we note that with mf, — 4.7 GeV and the cutoff Eq — 1.8 GeV 
employed in a recent analysis by the Belle Collaboration [6] one gets A = 1.1 GeV. For 
Eq = 1.6 GeV (a reference value adopted in [7, 27], which at present is below what can be 
achieved experimentally) one would obtain A = 1.5 GeV. As mentioned in the Introduction, 
in all previous analyses of the B —>■ Xgj decay rate it was assumed that, once Eq is taken 
below about 1.9 GeV, the sensitivity to hadronic physics essentially disappears, and the rate 
can be computed using a conventional OPE at the scale rrib. The main point of the present 
work is to show that this assumption cannot be justified, and that estimating theoretical 
uncertainties under the hypothesis that the expansion is in powers of ^^(mb) and Aqqb /n^b 
underestimates the magnitude of the true theoretical errors. As we will show, for values of Eq 
outside the shape-function region there are three relevant mass scales in the problem besides 
■'^QCD- They are the hard scale rrib, the hard-coUinear scale \/rr^K, and the low scale A 
itself. The values of these scales as a function of the photon-energy cutoff Eq are shown in 
Figure 1. The transition from the shape-function region to the region where a conventional 
OPE can be applied is not abrupt but proceeds via an intermediate region, in which a short- 
distance analysis based on a multi-stage OPE (MSOPE) can be performed. The transition 
from the shape-function region into the MSOPE region occurs when the scale A becomes 
numerically (but not parametrically) large compared with Aqcd- Then terms of order a" (A) 
and (Aqcd/ A)", which are non-perturbative in the shape-function region, gradually become 



Fcusp («.()")) In — + 7' («.()")) ■ 
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Figure 1: Dependence of the three scales fih = t^^h (sohd), /ij = \/rr^K 
(dashed), and iiq = ^ (dash-dotted) on the cutoff Eq, assuming m;, = 4.7 GeV. 
The gray area at the bottom shows the domain of non-perturbative physics. 
The hght gray band in the center indicates the region where the MSOPE should 
be applied. 

decent expansion parameters. Only for very low values of the cutoff {Eq < 1 GeV or so) it is 
justified to treat A and ^/nibA as scales of order rrib. 

Separating the contributions associated with these scales requires a multi-step procedure, 
which we develop in the present work. The first step, the separation of the hard scale from 
the intermediate scale, has already been achieved in (1). To proceed further we use two 
crucial recent developments. First, integrals of smooth weight functions F(u}) with the shape 
function S{u},/j) can be expanded in a series of forward i?-meson matrix elements of local 
HQET operators, provided that the integration domain is large compared with Aqcd [14, 15]. 
The reason is that the shape function can be written as the discontinuity of a two-point 
correlator in momentum space, and thus weighted integrals over S can be turned into contour 
integrals in the complex plane along a circle with radius set by the upper integration limit on 
uj (more precisely, a) — A). Specifically, the expansion takes the form [14] 

j^^'^dwSiuj^ti) Ficu) = 4^\A,i,) + K^/\A,i,) + . . . , (15) 

where K^^^ are calculable Wilson coefficient functions, A = uib — rrib and Ai are HQET 
parameters (which for the time being are defined in the pole scheme) [20], and the dots 
represent terms of order (Aqcd/^)^ or higher. Note that with A = m;, — 2i?o as defined above 
we have A + A = Mb — 2Eq = Ae, which coincides with the upper limit for the integration 
over a) in (1). The perturbative expansions of the coefficient functions K^^^ can be trusted as 
long as /X ~ A. In order to complete the scale separation, it is therefore necessary to evolve the 
shape function in (1) from the intermediate scale Hi ~ v^m^A down to a scale //q ~ A. This 
can be achieved using the analytic solution to the integro-differential RG evolution equation 
for the shape function in momentum space obtained in [14, 18]. These manipulations will be 
discussed in detail in the following two sections. 
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As a final comment, we stress that the main purpose of performing the scale separation 
using the MSOPE is not that this allows us to resum Sudakov logarithms by solving RG 
equations. Indeed, the "large logarithm" In(m5/A) ^ 1.5 is only parametrically large, but 
not numerically. What is really important is to disentangle the physics at the low scale 
/iq ~ A, which is "barely perturbative" , from the physics associated with higher scales, where 
a short-distance treatment is on much safer grounds. It would be wrong to pretend that all 
perturbative effects in B Xg'-f decays are associated with the short-distance scale mi, ^ 
^QCD- The MSOPE allows us to distinguish between the three coupling constants as{mh) ^ 
0.22, as{\/mbA) ^ 0.29, and ^^(A) pa 0.44 (for A = 1.1 GeV), which are rather different 
despite the fact that there are no numerically large logarithms in the problem. Given the 
values of these couplings, we expect that scale separation between A and nib is as important 
as that between nib and the weak scale Mw- 



3 Evolution of the shape function 



The renormalized shape function obeys the integro-differential RG evolution equation 

d 



din II 



dill' js{i^,i^',IJ') S{lj',h) , 



(16) 



where the anomalous dimension can be written in the form 



-fs{L0,L0',ll) = -2rcusp(Qis(A*)) 



W — LO' 



(17) 



This form was found in two recent one-loop calculations of the ultra-violet poles of non-local 
HQET operators [14, 15] . A brief history of previous investigations of the anomalous-dimension 
kernel can be found in the first reference. The structure of (17) was derived first by Grozin and 
Korchemsky [38] , who also computed the anomalous dimension and argued that the functional 
form of 7s(cj, cj', /i) shown above holds to all orders in perturbation theory. A conjecture for 
the two- loop expression of the anomalous dimension 7 is presented in Appendix A.l. 

The exact solution to (16) can be found using a technique developed in [18]. The equation 
is solved by the remarkably simple form 



S{u;,lii) = U2{iii,iio) 



-■yEV 



ii^q{u - (i^'y-'i 



where 



lnC/2(//i,/Xo) = 2S{iJ,o,fXi) + 2a^{fio,iii) , 77 = -2ar()Uo, /^i) ■ 



(18) 



(19) 



The functions S and a-p have been defined in (12). Similarly, the function is defined in com- 
plete analogy with ay, but with 7' replaced with the anomalous dimension 7 in (17). Explicit 
equations for these functions are given in Appendix A.l. The next-to-leading logarithmic ap- 
proximation to (18) was first derived in [14]. We note that a similar (but not identical) result 
was found in [39] based on a one-loop calculation of the anomalous-dimension kernel. 
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Relation (18) accomplishes the evolution of the shape function from the intermediate scale 
down to the low scale fiQ ~ A. When this result is inserted into the factorization formula (1), 
it is possible to perform the integrations over and uj analytically, leaving the integration 
over a)' until the end. Using the expression for the jet function in (9), we find that the leading 
contribution to the decay rate is given by 

(20) 

where 



Jo \A*o/[ 47r j ^ Mb — CO'' 

(21) 

The function ^3(77, S) is a special case of the polynomial 



The next-to-leading order corrections from the jet function are encoded in the operator 
J( A) = 2 (^In ^ + - [^h{ri) + 3] ^ + 

+ 2h\ri) + 3/1(77) - 2/i'(77) + 7 - ^ , (23) 

where 

/i(r/) =^(r/)+7E + - =^(l + r/)+7E (24) 

is the harmonic function generalized to non-integer argument, and the derivatives 5^ = d/dr] 
in (23) act on the function ^3(77, 5) in (21). Note that this result has a smooth limit for 
A*o IJ'i, in which case r] ^ 0, f/2(A*i,A*o) 1, h{r]) —>■ 0, /^'(r^) — >• 7r^/6, and we obtain an 
expression equivalent to the original result in (1). 



4 Short-distance expansion of the convolution integral 

The remaining task is to expand the integral over the shape function in (21) using an OPE, 
relating it to forward B-meson matrix elements of local HQET operators, as indicated in (15). 
As explained in [14], this can be done whenever A = A^; — A = rrih — 2Eq is large compared 
with Aqcd- For a given weight function F, the matching coefficients K^^^ are determined 
in the usual way by computing the integral in perturbation theory, expanding in powers 
of external momenta, and writing the answer as a linear combination of Wilson coefficients 
multiplying the matrix elements of local HQET operators. This matching calculation can be 
done using free partons in the external states and employing any infra-red regulator scheme 
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that is convenient. We use on-shell external heavy-quark states with residual momentum k 
chosen such that v-k = In this case, the perturbative expression for the renormalized shape 
function at one-loop order is [14, 15] 



parton('^, IIq) = 5{uJ - K + Tl ■ k) 1 



(//(-)) -K^ 
24 



TT 



TT 



1 



Cj — h. -\- n ■ k 



uj-K + n-k 
2 In h 1 



[mo] 



(25) 



Using this result one can perform the integral over the shape function in (21) analytically. 
The answer is then Taylor-expanded in powers of n • A;. The terms up to second order in 
this expansion are identified with the forward 5-meson matrix elements of the operators hh, 
h in ■ Dh, and h {in ■ Dyh, respectively, where n'^ = (1, 0, 0, 1) is a light-hke vector. The values 
of these matrix elements are given by 1, 0, and — Ai/3 [8]. They do not receive radiative 
corrections in the regularization scheme adopted here. Operators of dimension six or higher 
would mix under renormalization. Also, in order to find their Wilson coefficients it would be 
necessary to perform matching calculations with external gluon states. However, it will be 
sufficient for all practical purposes to truncate the expansion after the second term, keeping 
only operators of dimension up to five. The result of this calculation is 



I{Eo) = mH - 




^ ^ CV«^ ^^^^ ^ C,as{^^o) ^^^^ 



X 



PsiV, — ) + 



47r 

r]{r]-l) (-Ai) 
2 3A2 



47r 



+ ... 



(26) 



where J (A.) has been defined in (23), and 



5(A) = -4 fin — + +4:[2h{ri)-l](\n— + dr]-4:h\ri)+ih{ri)+ih'{ri)-^. (27) 

\ Ho J V /^O / D 

Wc have restricted ourselves to include only the leading power correction of order Ai/A^, 
dropping terms that are suppressed by additional powers of A/m^. This is necessary for 
consistency, because there exist other, unknown l/rrih and I /ml corrections from subleading 
shape functions, i.e., non-local HQET operators containing additional derivatives or insertions 
of soft gluon fields [40]. The Ai/A^ term is obtained by acting with (— Ai/6) 9^ on the leading- 
order term. According to (26), its effect can be included by simply adding a power correction 
to the function p:i{ri, d). 

The reader may ask why such an "enhanced" power correction was not found in previous 
analyses of the decay B — > Xg^, or of the related semileptonic decay B ^ Xulv. Common 
lore is that non-perturbative corrections to inclusive decay rates scale like {Aqcb / mb)^ and 
thus are very small. The reason is that so far power corrections in the OPE were computed 
at tree level only (an exception being [41]). While the terms displayed above have a non- 
zero leading-order coefficient after RG resummation, they vanish at tree level if the result is 



11 








-0.1 



0.6 



0.8 1 1.2 1.4 1.6 1.8 2 



A [GeV] 



Figure 2: Size of the enhanced power correction proportional to Ai/A^ in (26) 
relative to the leading term, as a function of A = m;, — 2Eo. 

expanded in fixed-order perturbation theory. Explicitly, we find to first order in ag 



tions had been computed beyond the tree approximation. 

Even though it is parametrically larger than the non-perturbative corrections from the 
conventional OPE in (2), the enhanced power correction in (26) remains small for all relevant 
values of A. This is illustrated in Figure 2, which shows the size of the power correction 
proportional to Ai to the function I{Eq) relative to the leading-order term. In the region 
of "perturbative" values of A, where the MSOPE can be trusted, the effect amounts to a 
reduction of the decay rate by less than 5%. It also follows that subleading power corrections 
of order Ai/(m6A) can safely be neglected. 

Equation (26), combined with (20), is our main result. Its numerical implications will be 
analyzed later, after including additional, small power-suppressed terms. A few comments are 
in order already at this point: 

1. The rate in (20) is formally independent of the three matching scales, at which we switch 
from QCD to SCET (/i/j), from SCET to HQET (yUj), and finally at which the non-local 
HQET matrix element (the shape-function integral) is expanded in a series of local 
operators (//o)- The explicit perturbative expressions for the functions H^{nh) in (4), 
J'(A) in (23), and 5(A) in (27) suggest that the "natural" choices for the three scales are 
Hh = fTT-i), yUj = y/rriiyA, and fiQ = A, as this removes all logarithms from these expressions. 
The latter two assignments are supported by the observation that, for a typical value 
T) 0.25, the coefficient function J7'(A) vanishes near /li = 1.08\/m(,A, while |5(A)| is 
minimized near //q = 1.16 A. Below, we will adopt the "natural" choices as our default 
values. In practice, a residual scale dependence arises because of the truncation of the 
perturbative expansion. Varying the three matching scales about their default values 
provides some information about unknown higher-order perturbative terms. 




(28) 



This effect would have shown up in the conventional heavy-quark expansion, if power correc- 
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2. In the limit where the intermediate and low matching scales //j and jdo are set equal 

to the hard matching scale fi^, our result reduces to the conventional formula used in 
previous analyses of the B — > Xg'j decay rate. However, this choice cannot be justified 
on physical grounds. 

3. After RG resummation the decay rate has a non-trivial dependence on the photon-energy 
cut Eq already at leading order in RG-improved perturbation theory and at leading power 
in A/mb, as reflected by the appearance of A^' in (26). 

4. In (20) we have accomplished a complete resummation of (parametrically) large loga- 
rithms at next-to-next-to-leading logarithmic order in RG-improved perturbation theory, 
which is necessary in order to calculate the decay rate with 0{as) accuracy. This is highly 
non-trivial in cases where Sudakov double logarithms are present. Specifically, it means 
that terms of the form a^L*^ with k = (n — 1), . . . , 2n and L — ln(mb/A) are correctly 
resummcd to all orders in perturbation theory. At a given order a", there arc (n + 2) 
such terms. To the best of our knowledge, a complete resummation at next-to-next-to- 
leading order has never been achieved before. For ease of comparison with the results of 
other authors, we provide in Appendix A. 2 an expansion of our result to second order 
in CKs, deriving the coefficients of the terms a'^L'^ with A; = 1, 2, 3, 4. 

5. Finally, we stress that the various next-to- leading order corrections in the expression for 
the decay rate obtained from (20) and (26) should be consistently expanded to order 
as before applying our results to phenomenology. Such next-to-leading order terms are 
contained in the functions Q^(/i/,), Wibdih), H^ifJ'h), Ui{iJ,h, fJ'i), t^2(Aii, A*o), V, and I{Eo). 
For instance, one should expand 



-"/ev p-ism 



1 _ (^£i _ «s(/^o) - (^silJ'i) ^^^^^ ^ 



(29) 



where 770 = ^ In ^"(^"j is the leading-order expression for 77 (see Appendix A.l). In 
practice, these expansions are readily automatized. 



5 Elimination of pole-scheme parameters 

The expression for the function I{Eq) in (26) has been derived under the implicit assumption 
that the 6-quark mass nib, the related parameter A = nib — and the HQET parameter Ai 
are defined in the on-shell scheme. While this is most convenient for performing calculations 
using heavy-quark expansions, it is well known that HQET parameters defined in the pole 
scheme suffer from infra-red renormalon ambiguities [42, 43, 44, 45]. As a result, the pertur- 
bative expansion in (26) would not be well behaved. It is thus necessary to replace the pole 
mass nib and the HQET parameter Ai in favor of some physical, short- distance parameters. 

For our purposes, the "shape-function scheme" defined in [14] provides for a particularly 
suitable definition of the heavy-quark mass and kinetic energy. A look at (21) shows that the 
pole mass actually never enters the expression for the decay rate. Rather, a factor {Mb — oj)^ 
appears under the integral over the shape function, which can be traced back to the factor 
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{Mb — P+y = SE^ in the original expression for the rate in (1). Roughly speaking, it is the 
average value of uj that determines the value of the difference {Mb — tj). This observation is 
the basis of the shape-function scheme. The idea is that a good estimate of the right-hand 
side of (21) can be obtained using the mean- value theorem, i.e., by replacing u with its mean 
value defined as 



J^'^diu S{lu,Ij,o) 



A(A, /io) = Mb - mb{A, //q) . 



(30) 



Here m5(A,/io) is the running shape-function mass defined in [14], which depends on a hard 
cutoff A in addition to the renormalization scale /^o- The quantity A in the shape-function 
scheme is defined by the implicit equation 



A = Ae - A(A, Ho) = mb{A, hq) - 2Eo . 



(31) 



(For simplicity, we write A instead of the more correct notation A(A, //o)-) Likewise, we define 
a kinetic-energy parameter /x^(A, /^o) via the variance of cu, 



{^'')a - {u:)\ 



S^''dwS{u),HQ) 



fJ'l{A,f^o) 



(32) 



The shape- function scheme provides a physical, short-distance definition of nib and /x^, which 
can be related to any other short-distance definition of these parameters using perturbation 
theory. The explicit form of these relations for some common renormalization schemes can be 
found in [14]. Here we need the relations to the parameters defined in the pole scheme. They 
are 



m 



pole 



mb{A,iio) + A 



TT 



1 - 2 In ■ 



I 2 /i^(A,/.o) A 
3 A2 



+ ... 



-Ai = /x^(A,//o) 



1 + 



CFOis{llQ) 



-3 In 



TT 



A 

1^0 



_^3^2 CFas{fJ'o) 



TT 



, A 

In \- ... . 

1^0 



(33) 



pole 



— 2 En. In order to 



The corresponding relation for Ap°^'^ follows from the fact that Ap°'^ = m 
introduce the parameters defined in the shape-function scheme, we perform these replacements 
in the expression for I{Eq) in (26) and expand the result consistently to order a^. (In the 
next-to-leading order terms we can simply replace the parameters of the pole scheme by the 
corresponding parameters of the shape-function scheme.) 

While the above choice appears most natural to us, it is by no means unique. For instance, 
we may avoid using the running quantities mb{lJL j , jj) and /j'^{/jf,fi) with "off-diagonal" scale 
choices fij ^ fihj using instead the parameters mi,{fi, fJ^) and ij,^{ij,, /j,), which are related to 
the parameters in the pole scheme by the simpler relations 



m 



pole 



, . CFOls{fJ') 

mb[iJ,,ii) +11 \- ... 



TT 



-Ai = i4{lJ',IJ') 



1 - 



CFas{lJ') 
27r 



+ 



(34) 



The parameter A is now determined by the equation A = mb{fJ^, fJ^) — 2Eq. The scale fi could 
naturally be taken to be /Iq. Alternatively, we may use the parameters of the shape-function 
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scheme defined at a fixed reference scale /x* = 1.5 GeV, at wliicli tlieir values have been 
determined to be mfe(/i„^,) = (4.65 ± 0.07) GeV and nKfi^, fi^) = (0.27 ± 0.07) GcV^ [14]. 
These determinations are based on various sources of phcnomcnological information, including 
T spectroscopy and moments of inclusive i?-meson decay spectra. In our numerical analysis 
in Section 8 we will present results for different variants of the shape-function scheme. 

6 Kinematic power corrections 

The results of the previous section provide a complete description of the B — > Xg^ de- 
cay rate at leading order in the l/mf, expansion, where the two-step matching procedure 
QCD —>■ SCET —>■ HQET is well understood. The matching coefficients and anomalous dimen- 
sions are known to the required order, so that the scale separation and RG resummation can 
be carried out with next-to-next-to-leading logarithmic accuracy. For practical applications, 
however, it is necessary to also include corrections arising at higher orders in the heavy-quark 
expansion. The leading non-perturbative corrections proportional to the HQET parameter Ai 
(or have already been included above. More important, however, are "kinematic" power 
corrections of order (A/m^)", which are not associated with new hadronic parameters. Un- 
like the non-perturbative corrections, these effects arise already at first order in A/m^, and 
they are numerically dominant in the region where A ^ Aqcd- Technically, the kinematic 
power corrections arise in the matching of QGD correlators onto higher-dimensional SGET 
and HQET operators. 

The corresponding terms are known in fixed-order perturbation theory, without scale sep- 
aration and RG resummation [26, 46] (see also [17]). To perform a complete RG analysis of 
even the first-order terms in A/m^ is beyond the scope of the present work. Since for typical 
values of Eq the power corrections only account for about 15% of the B Xg'-f decay rate, 
an approximate treatment will suffice. To motivate it, we note the following two facts: First, 
while the anomalous dimensions of the relevant subleading SCET and HQET operators are 
only known for a few cases [47], the leading Sudakov double logarithms are determined by 
the cusp anomalous dimension and thus are the same as for the terms of leading power. The 
reason is that the cusp anomalous dimension has a geometric origin. In the present case, 
it results from a product of time-like and light-like Wilson lines describing heavy and hard- 
collinear quark fields, respectively [48]. The leading Sudakov double logarithms are therefore 
the same as those resummed into the exponents S{iJ,h, l-ii) and 5'(//o, A*i) contained in the evo- 
lution functions Ui and U2 in (11) and (19). Secondly, all power-suppressed terms of order 
(A/m^)" are associated with gluon emission into the hadronic final state Xg. Because of the 
kinematic restriction to low-mass final states, i.e. Mx < Mb^e, the emitted gluon can only 
be hard-collinear or soft, but not hard. One should therefore associate a coupling as{Hi) or 
as{iio) with these terms. The leading power corrections are then of order as{ii)5\n5 with 
5 — A/mfc and ~ /Xj or /xq. After RG resummation, they can give rise to effects of order rjS, 
which are formally of zeroth order in the coupling constant. Not resolving the scale ambiguity 
for such terms introduces an uncertainty that is at most of order In^ 5. 

In order to at least partially account for resummation effects, we proceed as follows: We 
include the known power corrections from real gluon emission and associate the coupling as{iii) 
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with them. The Wilson coefficients Cj are evaluated at the hard scale jdh- We then multiply 
the answer with the evolution function Ui. This accounts correctly for the leading Sudakov 
logarithms in the evolution from the hard scale fih to the intermediate scale /ij. While the 
conventional parton-model calculation of the B — > Xg'y decay rate is performed with on-shell 
b quarks, we add a residual momentum such that pj, — rribV + k. In the light-cone component 
n ■ Pb — rrib + n ■ k we keep the n • k piece, because it is of the same order as the corresponding 
component n ■ p^,. of a hard-coUinear momentum. In all other components we neglect k. 
This accounts for some, but not all shape-function effects. The net result is that we must 
replace — > Mb — a) (and hence — > — a)) in the parton-model calculation, and then 
convolute the result with the leading-order shape function. In the approximation where the 
small parameter £ckm in (5) is set to zero (which is an excellent approximation given that we 
are dealing with power-suppressed effects), this yields 



Gla 



X 



C^a.(/..) ^ ^^.^^^^ ^^^^^^^ H^r^) - Re {C*M Cf^i^H 



i < j 



(35) 



The functions fij{S) vanish linearly with 6 and so are of order A/m^. They coincide with 
3/jj(5) in [17] except for the case of fyriS), which requires an additional subtraction due to 
the fact that the function P3{ri,6) in (26) already contains some power corrections resulting 
from the presence of the factor (Mb — P+)^ in (1). We find 



f77{S) = 3/77(<^) - 5 (12 In 5 + 9) + 52 ( 6 In (5 + y ) - 53 In (5 + ^) 



(36) 



The relevant expressions are 
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where xs — max(a;, 1 — 5), as previously z — {mc/m^'^, and 



G{t) 



-2arctan2^t/(4-t) 



; i <4, 



2 In 



(v^ + Vi^)/2j-y)" ; t>4. 



(38) 



The next step is to account for the evolution between yUj and /ig, and to evaluate the shape- 
function integrals for A Aqcd using the techniques described in Section 4. From (18) and 
(25), we find 



'-'parton('-t^, A'-ij — 'J2[IJ'i, fJ'O) ^/ x tm „ + 



(39) 



Tin) iil{u-KY-^ 

where A is defined in the shape-function scheme (see Section 5), and the dots represent terms 
of order Q;s(/io) and higher-order non-perturbative corrections, which we consistently neglect. 
Inserting this result into (35) yields 
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, (40) 



where 



dyrjy^-\l-ySfj\ 




yd 



(41) 



The definition of the "smeared" functions Fij(r],5) is such that Fij(0,5) — fij{S), Fij(r],0) — 
4(0) = 0, and F,,ir],l) = /,,{!). 

The result (40) has the desired features that the leading Sudakov double logarithms are 
correctly resummed in the product U1U2, and that the gluon-emission terms are associated 
with a low-scale coupling constant that is larger than as{nh)- However, we stress that while 
the result is correct when expanded in fixed-order perturbation theory to first order in a^, the 
resummation of single logarithmic terms is only approximate. After a complete RG resum- 
mation, terms of the form asln{A/mb), which arise from the In 6 terms in the expressions for 
the functions fij, would be resummed into functions of rj, e.g. 
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(42) 



The correct answer will contain more complicated functions of 77 as well as non-logarithmic 
next-to-leading-order corrections at the scales /ij and /io- 

While we expect that (40) gives a good approximation for the power-suppressed contri- 
butions to the B — > Xg'j decay rate, it would be important and conceptually interesting to 
explore the structure of power corrections further, using the effective field-theory technology 
developed here and in [14]. It should be possible (with a significant amount of work) to re- 
solve the scale ambiguity for the first-order power corrections in A /nib. Also, an effective 
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field-theory analysis would allow a more rigorous description of certain non-perturbative ef- 
fects, such as the A2/m^ term in (35), which models a long-distance contribution related to 
charm-penguin diagrams [49, 50], or the logarithmic mass singularity regularized by in the 
expression for fss in (37), which is related to fragmentation effects [51]. More generally, such 
an analysis would provide a transparent power counting for any long-distance contributions 
involving soft partons (not only heavy quarks) in the MSOPE. 



7 Ratios of decay rates 

The contributions from the three different short-distance scales entering our central result (20) 
and the associated theoretical uncertainties can be disentangled by taking ratios of decay rates. 
Some ratios probe truly short-distance physics (i.e., physics above the scale /ih ~ nT'b) ^-nd so 
remain unaffected by the new theoretical results obtained in this paper. For some other ratios, 
the short-distance physics associated with the hard scale cancels to a large extent, so that one 
probes physics at the intermediate and low scales, irrespective of the short-distance structure 
of the theory. These ratios are important, because they are insensitive to New Physics and 
just probe the interplay of hard-coUinear and low scales in the process. Below, we investigate 
examples of both classes of ratios. 





NP 


-f^7(/i/i) NP 




SM 





7.1 Ratios insensitive to low-scale physics 

Most importantly, physics beyond the Standard Model may affect the theoretical results for 
the B — > Xg^ branching ratio and CP asymmetry only via the Wilson coefficients of the 
various operators in the effective weak Hamiltonian. (An exception are unconventional New 
Physics scenarios with new light particles, such as a supersymmetric model with light gluinos 
and b squarks considered in [52].) As a result, the ratio of the B — Xgj decay rate in a 
New-Physics model relative to that in the Standard Model remains largely unaffected by the 
resummation effects studied in the present work. Prom (20), we obtain 

-I- power corrections. (43) 

The power corrections would introduce some mild dependence on the intermediate and low 
scales fii and fiQ, as well as on the cutoff Eq. 

Another important example is the direct CP asymmetry in i? — > Xg'y decays, for which we 
obtain 

Acp = '-—^ = m2^\Tj\ M2 + PO^®^ corrections, 44 

where H^{^h) is obtained by CP conjugation, which in the Standard Model amounts to re- 
placing £cKM ^CKM (4)- It follows that the predictions for the CP asymmetry in the 
Standard Model and various New Physics scenarios presented in [29] remain largely unaffected 
by our considerations. 
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7.2 Ratios sensitive to low-scale physics 

The multi-scale effects studied in this work result from the fact that in practice the B — > Xg^ 
decay rate is measured with a restrictive cut on the photon energy. As we have pointed out, 
this introduces sensitivity to the scales /li ~ ^/rnbA and /xq ~ A = mf, — 2Eo in addition to the 
hard scale fih ~ 'm-b- These complications would be absent if it were possible to measure the 
fully inclusive rate. It is convenient to define a function F{Eq) as the ratio of the B — 
decay rate with a cut divided by the total rate, 



F{E^) = 



(45) 



Because of a logarithmic soft-photon divergence for very low energy, it is conventional to define 
the "total" inclusive rate as the rate with a very low cutoff E^ = mf,/20 [17]. The denominator 
in the expression for F{Eq) can be evaluated using a conventional OPE, which corresponds 
to setting all three matching scales equal to Hh- The numerator is given by our expression in 
(20), supplemented by the power corrections in (40). We obtain 
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where S — A./mb, 5* = 1 — 2E^/mb — 0.9, and 
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+ 21n5* + 



2; "iblml ■ 



(47) 



The result (46) is RG invariant and so (formally) independent of the three matching scales /i/^, 
/^i, and /io, and at leading power it is insensitive to the hard matching corrections contained in 
H^{Hh). To an excellent approximation, the fraction function F{Eo) therefore applies to the 
Standard Model as well as to any New Physics scenario. Note also that the 6-quark mass enters 
the expression for the fraction function only at the level of power corrections. The prefactor 
mlrnKfih), which multiplies the total decay rate, cancels out in the ratio (45). Finally, we 
stress that the expression for F{Eq) given above still refers to the pole scheme. It is necessary 
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to eliminate the pole-scheme parameters mi, and Ai in favor of physical parameters before 
using this result. 

Another important example of a ratio that is largely insensitive to the hard matching 
contributions is the average photon energy defined as 



dE^ 



Eo dEy 

which has been proposed as a good way to measure the 6-quark mass or, equivalently, the 
HQET parameter A [53, 54]. The impact of shape-function effects on the theoretical prediction 
for this ratio has been studied in [17, 55] and was found to be significant. Here we study the 
average photon energy in the MSOPE region, where a model-independent prediction can be 
obtained. It is structurally different from the one obtained using the conventional OPE in the 
sense that contributions associated with different scales are disentangled from each other. We 
find (with 5 = A/ nib) 
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where P4{7], 6) is defined in (22), and 

d,ji5)^ fdxfij{x)-Sfij{S). (50) 



Analytical expressions for the functions dij{6) are given in Appendix A. 3. They vanish quadrat- 
ically for 5 — > and so give very small contributions for realistic values of the cutoff. We there- 
fore do not include RG resummation effects for these terms. The non-perturbative corrections 
involving the parameters Ai and A2 are taken from [4]. Note that the expression in brackets 
is a purely perturbative result free of hadronic parameters. When expanded in fixed-order 
perturbation theory, our result (49) reduces to an expression first obtained in [17]. 

We stress that the hard scale Hh ~ mj, affects the average photon energy only via second- 
order power corrections. This shows that it is not appropriate to compute the quantity (E^) 
using a simple heavy-quark expansion at the scale m^, which is however done in the con- 
ventional OPE approach [53, 54]. This observation is important, because information about 
moments of the B X^^ photon spectrum is sometimes used in global fits to determine the 
CKM matrix element iV^cfel along with HQET parameters. Keeping only the leading power 
corrections, which is a very good approximation, the above expression simplifies to 



2 2(1 + r;) 



,(l + ,) + ^£^fln!I^-M,)-3 1 



TT V A*0 2 l+T] 



+ .... (51) 
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Table 1: Compilation of input parameters entering the numerical analysis. 
The top-quark mass enters the expressions for the Wilson coefficients Cj. The 
strange-quark mass is required as an infra-red regulator in (37). Only the real 
part of £cKM is needed for the calculations in this work. 



Jr^arameter 


Value 


Source 


mb{fi^,IJ^) [GcV] 


4.65 ±0.07 


[14] 


[GeV^] 


0.27 ±0.07 


[14] 


mb{mb) [GeV] 


4.25 ±0.08 


[56] 


fnc(rnc) [GeV] 


I. ZD ± O.iO 


[57] 


mr'" [GeV] 


178.0 ±4.3 


[58] 


rris/mh 


0.02 


[17] 


[ps] 


1.604 ±0.016 


[57] 




0.1187 ±0.0020 


[57] 




40.4 tJ:^ 


[59] 


Re(£cKM) [10-3] 


9.8^:1 


[59] 


Ai [GeV^] 


-0.25 ±0.20 


[60, 61] 


A2 [GeV2] 


0.12 





In this approximation, (E^) only depends on physics at the intermediate and low scales Hi and 
IIq. The next-to-leading order perturbative corrections in this formula are numerically quite 
significant. For £"0 = 1.8 GeV, and taking the default scale choices fii = \/miA and /iq = A, we 
find in the pole scheme (using m^°'^ = 4.8 GeV for the purpose of illustration) (E^) ^ [2.27 + 
0.29as{y/rnbA) — 0.19^5 (A)] GeV. Eliminating the pole mass nib in favor of the 6-quark mass 
mb{A., A) defined in the shape-function scheme, we obtain {E^) [2.222 + 0.254Q;s(\/m;,A) + 
0.009q;s(A)] GeV ^ 2.30 GeV. When the 6-quark mass is defined in the shape-function scheme, 
the average photon energy is numerically very close to |mfe(A, A) Ri 2.33 GeV, meaning that 
the first term in (51) dominates. Note also that the correction proportional to the low-scale 
coupling q;s(A) is largely reduced in this scheme, ensuring an improved perturbative behavior. 

8 Numerical results 

We are now ready to present the phenomenological implications of our findings. Table 1 
contains a list of the input parameters entering the analysis together with their present uncer- 
tainties. We have inflated the error on Ai obtained by averaging the values quoted in [60, 61] 
from 0.06 GeV^ to 0.20 GeV^, taking into account that this parameter is affected by infra-red 
renormalon ambiguities [44, 45]. We vary the quark masses mg and irib independently, in which 
case ^/z — mc{iJih) / f^bil^h) — 0.221 ± 0.027. Additional uncertainties related to the possibility 
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that the proper normahzation point for the charm-quark mass in penguin loop graphs may be 
significantly lower than the hard scale fih arc considered part of the perturbative error. 

The most important correlations between input parameters arc implemented as follows. 
We consider the two 6-quark masses m(,(/x*, /i.*) and m(,(mft) as being fully correlated and vary 
their values simultaneously. The same applies to the values of the parameters and 
Ai. Next, we use that the value of mi, is strongly anti-correlated with that of \Vcb\, because 
the most precise determination of is obtained from the analysis oi B X^liy decay 
distributions. A recent study in [62] quotes a correlation coefficient c = —0.49 between rrif, 
and \Vcb\- CKM unitarity ensures that \Vcb\ is to a very good approximation equal to the 
product |V^*V^{,|, so that the same anti-correlation can be assumed between nib and iV^^^bl- 

Before presenting our results, we reiterate that to apply the formulae derived in this work we 
must first eliminate the parameters rub and Ai defined in the pole scheme in terms of physical 
parameters defined in the shape-function scheme and then expand the answer consistently 
to 0{as), treating ratios such as as{ni) / as{nh) and as{no) / as{ni) as 0{1) parameters. This 
expansion is readily automatized. Throughout, we use the 3-loop expression for the running 
couphng asi/ji) defined in the MS scheme [57]. 



8.1 Partial B — > Xg'y branching ratio 

We begin by presenting predictions for the CP-averaged B — > Xg'y branching fraction with a 
cutoff > Eq applied on the photon energy measured in the S-meson rest frame. Lowering 
Eq below 2GeV is challenging experimentally. The first measurement with Eq — 1.8 GeV has 
recently been reported by the Belle Collaboration [6]. It yields^ 

Br(S ^ X,7) = (3.38 ± 0.30 ± 0.29) • 10"^ , 

^ * " £;o=1.8GeV ^ ' ' 

{E^) = (2.292 ± 0.026 ± 0.034) GeV . (52) 

So— 1-8 GeV 

For £"0 = 1.8 GeV we have A 1.1 GeV, which is sufficiently large to apply the formalism 
developed in the present work. We will also present results for Eq — 1.6 GeV because this 
value has been used in some theoretical studies, although it has not yet been achieved in 
an experiment. (For comparison, the value Eq = 2.0 GeV adopted in the CLEO analysis [5] 
implies A ^ 0.7 GeV, which we believe may be too low for a short-distance treatment.) 

We first set all input parameters to their default values and study the dependence of 
the branching ratio on the three matching scales Hh, l-ii, and /iq. The sensitivity of our 
predictions to variations of the matching scales provides an estimate of unknown higher-order 
perturbative corrections. We shall study three different version of the shape-function scheme 
for the definition of the 6-quark mass and the kinetic-energy parameter /x^, as discussed in 
Section 5. In the first scheme (called "RS 1") we use the parameters mb(A,/io) and /i^(A,/io) 
defined in (33). In the second scheme ("RS 2") we instead use mb{iiQ, /iq) and /x^(/xo,A*o) from 
(34). Finally, in the third scheme ("RS 3") we employ the parameters mb{[x*, n*) and //^(/x*, /i*) 
renormalized at a fixed scale //* = 1.5 GeV, at which their values have been determined in 

^To obtain the first result we had to undo a theoretical correction accounting for the effects of the cut 
> 1.8 GeV, which had been applied to the experimental data. 
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Table 2: B ^ Xg^y branching ratio with estimates of perturbative uncertainties 
obtained by variation of the matching scales, for three variants of the shape- 
function scheme. See text for explanation. 



Eq 



Scheme 



Br [10- 



lii jJLQ Sum 



Power Cors. 



Combined 



1.8 GeV 



RS 1 
RS 2 
RS 3 



3.37 
3.38 
3.36 



+0.02 
-0.00 
+0.02 
-0.00 
+0.02 
-0.00 



+0.25 
-0.37 
+0.25 
-0.37 
+0.25 
-0.37 



+0.41 +0.48 
-0.03 -0.37 



h0.15 
-0.18 
h0.18 
-0.18 



hO.29 
-0.41 
h0.30 
-0.41 



+0.12 

-0.07 
+0.12 
-0.07 
+0.12 
-0.07 



+0.49 
-0.38 
+0.31 
-0.42 
+0.32 
-0.42 



1.6 GeV 



RS 1 
RS 2 
RS 3 



3.47 
3.47 
3.48 



+0.02 
-0.00 
+0.02 
-0.00 
+0.02 
-0.00 



+0.28 
-0.39 

+0.28 
-0.39 
+0.28 
-0.39 



+0.14 
-0.01 

+0.13 
-0.14 

+0.18 
-0.13 



+0.31 
-0.39 

+0.31 
-0.41 

+0.33 
-0.41 



+0.10 
-0.05 

+0.10 
-0.05 
+0.10 
-0.05 



+0.33 
-0.39 

+0.33 
-0.41 

+0.34 
-0.41 



[14]. In the schemes RS 1 and RS 2, these reference values are evolved to other scales using 
equations derived in [14]. 

The matching scales are independently varied about their default values /ih — iTib, Hi — 
^^uibA, and Hq — A by multiplying them with factors between 2/3 and 3/2. Thus, for 
nib = 4.7GcV and Eq = 1.8 GeV, we vary e [3.13, 7.05] GeV, /x^ £ [1.52, 3.41] GcV, and 
/io G [0.73, 1.65] GeV, while for Eq = 1.6 GcV the latter two ranges are replaced by /ij G 
[1.77, 3.98] GeV, and /Iq e [1.0, 2.25] GeV. Together, this covers a conservative range of scales. 
The resulting variations of the branching ratio are shown in Table 2. 

We observe an excellent stability of our predictions with respect to variations of the hard 
matching scale fih- In fact, the sensitivity is so small that it cannot reasonably be taken as an 
indication of the size of higher-order terms in the expansion in powers of as{lJLh)- The sensitivity 
to variations of the intermediate matching scale Hi is more pronounced. The numbers suggest 
that terms of order a1{iii) could impact the branching ratio at the 10% level, which appears 
entirely reasonable given that Q;s(/Xj) ~ 0.3. The sensitivity to the low matching scale //q turns 
out to be rather small. The coefficient of the ^^(/io) term depends on the scheme adopted 
for the definition of the parameters rrib and /x^, and it appears that in the three schemes 
considered here this coefficient is numerically small. While it is not guaranteed that this 
feature will persist in higher orders, the observation of good stability at the scale /xq suggests 
that the shape-function scheme captures the most important low-scale effects and absorbs 
them into the running &-quark mass and the parameter /x^. The column labeled "Sum" shows 
the combined uncertainty obtained by adding the three scale variations in quadrature. The 
next column, labeled "Power Cors." , gives an estimate of the perturbative uncertainty in our 
treatment of kinematic power corrections, as discussed in Section 6. It is obtained by studying 
two variants of the expression (40), one where we set ps — > 1 and Fij — > /jj, and one where in 
addition we neglect all anomalous dimension functions except those governed by Pcusp- In both 
cases, we obtain expressions that differ from (40) by terms that are beyond the accuracy of 
our calculation. The resulting changes in the branching ratio arc the same in all schemes and 
range between 1.5 and 3.5%, corresponding to a 10-25% uncertainty in the size of the power- 
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Table 3: B ^ Xg'y branching ratio with estimates of theoretical uncertainties 
due to input parameter variations as listed in Table 1. The upper (lower) sign 
refers to increasing (decreasing) a given input parameter. 





Br [10-^] 


rub rUc mt V^yth tb 0(s{Mz) 


Combined 


1.8 GeV 


3.38 


trsl ±0.10 ±0.04 tut ±0.03 tf^l 


+0.32 
-0.30 


1.6 GeV 


3.48 


±0.04 tlf, ±0.03 ±0.10 


+0.32 
-0.29 



suppressed contributions themselves. Finally, the last column in the table shows our estimates 
for the total perturbative uncertainty in the prediction of the branching ratio, which we find 
to be of order 10%, significantly larger than previous estimates. For example, the authors of 
[7, 27] argued in favor of a total perturbative error of only 4%. 

The remaining uncertainties in our predictions are due to input parameter variations. They 
are essentially the same in the three renormalization schemes and are summarized in Table 3 
for the case of RS 2. The last column shows the combined errors, added in quadrature. 
They are dominated by the uncertainties in the 6-quark mass and in \Vts\, whose significant 
anti-correlation (c = —0.49) is taken into account in computing the total error. Parameter 
dependences not shown in the table have a negligible effect (< 1%) on the branching ratio. 
Note that in contrast to previous authors we do not divide the theoretical expression for 
the B Xg'-f decay rate by a semileptonic rate, but present an absolute prediction for the 
branching ratio itself. Once the correlation between parameters is properly taken into account, 
normalizing T{B — > Xg'y) to the semileptonic rate r{B Xlv) does not lead to a significant 
reduction of the theoretical uncertainties. 

The above results can be combined into the new Standard Model predictions 

Br(5 ^ ^.7)|^^^,_3e^^ = (3.38 [pert.] t'^fo [pars.]) x 

Br(S ^ ^.7)|^^^^ = (3.471^:1? [pert.] t',^ [pars.]) x 10"^ , (53) 

where we use the mass renormalization scheme RS 2 as our default. The first error refers 
to the perturbative uncertainty and the second one to parameter variations. The first value 
is in excellent agreement with the experimental result (52). Comparing the two results, and 
naively assuming Gaussian errors, we find that^ 

Br(5 ^ X37)exp - MB ^ ^s7)sM < 1-3 ■ 10"^ (95% CL) . (54) 

We stress that, mainly as a result of the enlarged theoretical uncertainty but also due to the 
use of more recent data, this bound is much weaker than the one derived in [7]. where this 
difference was found to be less than 0.5 ■ 10^'^. Consequently, we obtain weaker constraints on 



^We do not use the CLEO data [5] in deriving this bound, because the choice Eq = 2 GeV does not allow 
for a model- independent treatment of the effects of the cut. 
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Table 4: S — > event fraction F{Eq) with estimates of perturbative uncer- 
tainties obtained by variation of the matching scales, for three variants of the 
shape-function scheme. See text for explanation. 



Eq 


Scheme 


F{Eo) [% 








Ato 


Sum 


Power Cors. 


Combined 


1.8 GeV 


RS 1 


89.1 


+2.4 
-2.2 


+1.4 
-5.0 


+2.4 
-5.4 


+3.6 
-7.7 


+4.7 
-2.4 


+5.9 
-8.1 




RS 2 


89.1 


+2.4 
-2.3 


+1.5 
-5.0 


+2.6 

-4.0 


+3.8 
-6.8 


+4.7 
-2.4 


+6.0 
-7.2 




RS 3 


89.2 


+2.5 
-2.3 


+1.3 
-5.0 


+2.5 
-3.7 


+3.8 
-6.6 


+4.6 
-2.4 


+6.0 
-7.0 


1.6 GeV 


RS 1 


93.1 


+2.8 
-2.6 


+2.7 
-5.7 


+2.6 
-2.3 


+4.7 
-6.7 


+3.9 
-1.9 


+6.1 
-7.0 




RS 2 


93.1 


+2.8 
-2.6 


+2.7 
-5.7 


+2.4 
-2.5 


+4.6 
-6.8 


+3.9 
-1.9 


+6.0 
-7.1 




RS 3 


93.1 


+2.8 
-2.6 


+2.7 
-5.7 


+2.6 
-2.3 


+4.7 
-6.7 


+3.9 
-1.9 


+6.1 
-7.0 



New Physics parameters. For instance, for the case of the type-II two-Higgs-doublet model, 
we may use the analysis of [63] to obtain the bound 

mH+ > approx. 200 GeV (95% CL) , (55) 

which is significantly weaker than the constraints m//+ > 500 GeV (at 95% CL) and mH+ > 
350 GeV (at 99% CL) found in [7]. To find the precise numerical value for the bound would 
require a dedicated analysis, which is beyond the scope of this paper. 

8.2 Event fraction ^(^^o) 

As an alternative way to discuss the effects of imposing the cut on the photon energy, we 
study the fraction function F{Eo) defined in (45), which up to power corrections is insensitive 
to the short-distance physics encoded in the Wilson coefficients Cj. The sensitivity of F{Eq) 
to scale variations is studied in Table 4, which is analogous to Table 2 for the branching ratio. 
We find that the fraction function exhibits a stronger sensitivity to the hard scale iih than 
the branching ratio, changing by about 3% as is varied between Im^j?) and ?>mi,/2. The 
sensitivity to variations of the matching scales //^ and //q follows the same pattern as in the case 
of the branching ratio, but the variations arc somewhat smaller in magnitude. Note that there 
is a difference between the function F{Eq) and the branching ratio as far as the dependence on 
fiQ is concerned, because the factor ml present in (26) and (40) cancels in the ratio (46). Since 
in the shape-function scheme the pole mass is expanded in a series in ^^(/xo), this has an 
effect on the perturbative expansion. Finally, the perturbative uncertainties in the calculation 
of the power-suppressed terms are again at the level of a few percent. Our estimate for the 
combined perturbative error is presented in the last column. 

In contrast to the B branching ratio, the fraction function F{Eq) is independent of 

several input parameters (i.e., mb{mb), \VflVtb\, tb, \i,2, and £ckm), and it shows a very weak 
sensitivity to variations of the remaining parameters. This is illustrated in Table 5, which 
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Table 5: S — > Xg'j event fraction F{Eo) with estimates of theoretical uncer- 
tainties due to input parameter variations as listed in Table 1. The upper 
(lower) sign refers to increasing (decreasing) a given input parameter. 







nil, III, (\,{Mz) III 


Coml)iiicd 


1.8 GeV 


89.1 


1?:^ ±0.5 to.3 


+1.3 
-1.6 


1.6 GeV 


93.1 


±0.3 ±0.4 ^0.5 ^0.1 


+0.7 
-0.8 



summarizes the resulting theoretical uncertainties for the case of RS 2. The combined errors 
are of order 1% and thus almost negligible compared with the perturbative uncertainties. 
In summary, we obtain 

F(1.8GeV) = (89 1^ [pert.] ± 1 [pars.])% , 

F(1.6 GeV) = (931? [pert.] ± 1 [pars.])% . (56) 

This is the first time that these fractions have been computed in a model independent way. 

The result corresponding to Eq = 1.8 GeV may be compared with the values {95.Sll2^)% and 
(95 ± 1)% obtained from the study of shape-function models in [17] and [55], respectively. In 
these studies, perturbative uncertainties have been ignored. A calculation in the conventional 
OPE approach gives a similar result, (95.2l2!9)% [7], where the authors took the error estimate 
from [17] . In the present work, we obtain a smaller central value with a larger uncertainty. 

As mentioned in the Introduction, the fraction function F[Eo) can be used to combine 
our study of multi-scale effects with other, independent calculations of the total B Xg'j 
branching ratio, both in the Standard Model and in extensions of it. For instance, we may use 
the result (3.70 ± 0.31) x 10~^ for the total branching ratio in the Standard Model obtained 
from [7, 27] (where the error contains a 4% perturbative uncertainty) and combine it with 
(56) to find 

Br{B ^ ^.7)|^^^,_3Q^v = (3.30 [pert.] ± 0.28 [pars.]) x 10"^ , 

Br{B ^ ^.7)|^^^^ = (3.441°:^,^ [pert.] ± 0.27 [pars.]) x 10"^ . (57) 

Compared with (53), the perturbative uncertainty is reduced by about 15-25%. However, only 
the reduction in the /iq dependence can be taken seriously, as the /li dependence is formally 
the same in (53) and (57). The insignificant reduction of the parameter uncertainties is partly 
due to the fact that the authors of [7, 27] take smaller parameter variations than those in 
Table 1, namely ±100 MeV for rric, ±30MeV for mh, and an error on the ratio | V^^Vj;,!/! V^&| 
that is half as big as what one would obtain using the global fit results compiled in [59] . 

As a final remark, we compare our results for the branching ratio with a cut at = 1-6 GeV 
in (53) and (57) with the benchmark value (3.57 ± 0.30) x 10~^ corresponding to the most 
recent calculation [27] published prior to the present work. Our central values are about 3% 
lower and, more importantly, the theoretical uncertainties we find are about 50% larger. 
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Table 6: Scale dependence and parameter variations for the average photon 
energy in B ^ Xg^ decays. See text for explanation. 



{£,) [MeV] 


/'/, /'/: /'o 


Conil). 


/"(;(/'*■/'*) -Ai (\,{Mz) 


Conil). 


2272 


1-1 +19 +48 
='=-'- -17 -70 


+51 
-72 


±37 ±10 +1 


±39 



8.3 Average photon energy 

The last quantity we wish to explore is the average photon energy. As discussed in Section 7.2, 
this quantity is almost insensitive to high-scale physics as well as to non-perturbative hadronic 
effects. However, it is very sensitive to the interplay of physics at the intermediate and low 
scales, as illustrated by the approximate relation (51). Our predictions for (E^) and its 
theoretical uncertainties are summarized in Table 6 for the case Eq = 1.8 GeV, corresponding 
to the cut employed in [6]. Since in this case the differences between the three variants of 
the shape-function scheme are insignificant, we only show results for RS 2. As expected, 
we find essentially no dependence on the hard matching scale, a modest dependence on the 
intermediate scale, and a more pronounced sensitivity to the low scale. The combined errors 
from scale variations are of order 50-70 MeV. The study of uncertainties due to parameter 
variations exhibits that the prime sensitivity is to the 6-quark mass, which is expected, since 
(E^) = m{,/2 + . . . to leading order. The next- important contribution to the error comes from 
the HQET parameter Ai. The total error is about 40MeV. 

Combining these results, we have to a very good approximation 

(^^)U...o.v = P.2T;S:S?)GeV+^-^, (58) 

where the error accounts for the perturbative uncertainty. The central values for the relevant 
input parameters are m{,(//*,/x*) = 4.65 GeV and Ai = — 0.25GeV^, and the quantities Srrib 
and SXi parameterize possible deviations from these values. Our prediction is in excellent 
agreement with the Belle result in (52). This finding provides support to the value of the 
6-quark mass in the shape- function scheme extracted in [14]. We stress, however, that the 
large perturbative uncertainties in the formula for (E^) impose significant limitations on the 
precision with which rub can be extracted from a measurement of the average photon energy. 
Our estimate above implies a perturbative uncertainty of 5m{,[pcrt.] = 1 Jqo McV in the ex- 
tracted value of rrib, which could only be reduced by means of higher-order calculations. This 
uncertainty is in addition to twice the experimental error in the measurement of {E^), which 
at present yields 6mi,[exp.] ~ 86MeV. 
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9 Conclusions and outlook 



In this work, we have performed the first systematic analysis of the inclusive decays B — > Xg'y 
in the presence of a photon-energy cut > Eq, where Eq is such that A — mj, — 2Eq can 
be considered large compared to Aqcd, while still A <^ mfe. This is the region of interest to 
experiments at the B factories. The first condition (A ^ Aqcd) ensures that a theoretical 
treatment without shape functions can be applied. However, the second condition (A ^ m;,) 
means that this treatment is not a conventional heavy-quark expansion in powers of asirub) 
and Aqcd /rrib. Instead, we have shown that three distinct short-distance scales are relevant 
in this region. They are the hard scale nib, the hard-coUinear scale y/iribA., and the low scale 
A. To separate the contributions associated with these scales requires a multi-scale operator 
product expansion (MSOPE), wc which have constructed in this work. 

Our approach allows us to study analytically the transition from the shape-function re- 
gion, where A ~ Aqcd, into the MSOPE region, where Aqcd -C A ^ m^, into the region 
A = 0{mb), where a conventional heavy-quark expansion applies. This is a significant im- 
provement over previous work. For instance, it has sometimes been argued that exactly where 
the transition to a conventional heavy-quark expansion occurs is an empirical question, which 
cannot be answered theoretically. Our formahsm provides a precise, quantitative answer to 
this question. In particular, for B Xsj with realistic cuts on the photon energy one is 
not in a region where a simple short-distance expansion at the scale nib can be justified. The 
precision that can be achieved in the prediction of the B — > Xg^ branching ratio is, ultimately, 
determined by how well perturbative and non-perturbative corrections can be controlled at 
the lowest relevant scale A, which in practice is of order 1 GeV. Consequently, we find larger 
theoretical uncertainties than previous authors. These uncertainties are dominated by yet 
unknown higher-order perturbative effects. Non-perturbative, hadronic effects at the scale A 
appear to be small and under control. 

Our treatment of the B Xs'~f branching ratio includes a complete resummation of 
logarithms ln(A/mb) at next-to-next-to- leading order in renormalization-group improved per- 
turbation theory. This level of precision has not been achieved before. Besides the calcu- 
lations performed here and in [14, 18], we have used multi-loop calculations for the cusp 
anomalous dimension [33, 34], the anomalous dimension of the shape function [64] (which 
we have corrected, see Appendix A.l and also [65]), and the anomalous dimension of the 
leading-order current operator in soft-collinear effective theory [37]. These ingredients are 
needed in order to achieve a complete separation of the perturbative corrections controlled 
by the three couplings as{nib), Q;s(v^mbA), and as{A), which differ in magnitude by about 
a factor 2. Our prediction for the CP-averaged B Xg'y with a cut Eq = 1.8 GeV is 
Bt{B X,7) = (S.SStoil [pert.] toU [pars.]) x lO"! With this cut (89 ± 1)% of all events 
are contained. The theory uncertainty wc estimate is significantly larger than that found by 
previous authors, and this fact has important implications for searches of New Physics in 
radiative B decays. Quite generally, the constraints on model parameter space have to be 
relaxed significantly. We have illustrated this fact with the example of the typc-II two-Higgs- 
doublet model, for which we find that the lower bound on the charged-Higgs mass is reduced 
to approximately 200 GeV. 

This is not the first time in the history of B ^ Xsj calculations that issues of scale setting 
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have changed the prediction and error estimate for the branching ratio. In [22], Czarnecki and 
Marciano have pointed out that the electromagnetic coupHng a in the expression for the decay 
rate should be identified with the fine-structure constant (normalized at = 0), and not with 
a{m^ renormalized at the scale of the heavy quark in the decay. This lowered the prediction 
for the branching ratio by about 5%. More recently, Gambino and Misiak have argued that the 
charm-quark mass, which enters the next-to-leading order corrections to the B Xg^ rate 
via penguin loops, should be identified with a running mass mc(/u) with /i ^ mh rather than 
with the pole mass [7]. This observation increased the prediction for the branching ratio by 
about 8%, and at the same time it increased the error estimate associated with the value of the 
ratio rric/mb, which before had been taken to be the (rather well known) ratio of the two pole 
masses. The point we emphasize in the present work, namely that some effects in S — > Xs^ 
decays should be described by the couplings as{VrnbA) and q;s(A) (and power corrections at 
the scale A) rather than as{mb), is of a similar nature. However, in our case the change in 
perspective about the theory of 5 — > Xgj decay is more profound, as it imposes limitations 
on the very validity of a short-distance treatment. If the short-distance expansion at the 
scale A fails, then the rate cannot be calculated without resource to non-perturbative shape 
functions, which would introduce an irreducible amount of model dependence. In practice, 
while A ^ 1.1 GeV (for Eq ^ 1.8 GeV) is probably sufficiently large to trust a short-distance 
analysis, it would be unreasonable to expect that yet unknown higher-order effects should be 
less important than in the case of other low-scale applications of QCD, such as in hadronic r 
decays. 

Given the prominent role of B ^ Xg^ decay in searching for physics beyond the Standard 
Model, it is of great importance to have a precise prediction for its rate in the Standard Model. 
The present work shows that the ongoing effort to calculate the dominant parts of the next-to- 
next-to-leading corrections in the conventional heavy-quark expansion is only part of what is 
needed to achieve this goal. Equally important will be to compute the dominant higher-order 
corrections of order ^^(A) and Q;g(-\/mbA), and to perform a renormalization-group analysis of 
the leading kinematic power corrections of order A/m^. In fact, our error analysis suggests that 
these effects are potentially more important that the hard matching corrections at the scale 
mij. Let us finish by mentioning two possible approaches for addressing the issue of higher- 
order perturbative effects at the intermediate and low scales: First, it would be interesting to 
calculate the terms of order PqoI at the scales /i, x/m^A and /ig ~ A. While this would fall 
short of a complete calculation of 0{a1) corrections, the "BLM terms" associated with the j3 
function are often numerically dominant [66, 67]. We stress that the known C(/3oQ;^) terms 
computed in the conventional heavy-quark expansion [54] are not sufficient for this purpose. 
Separate computations of 0(/3oQ;f) terms at the scales ~ rubi A*i ~ V^^fcA, and //q ~ A 
would be required to perform a meaningful BLM scale setting. This statement is explained 
in Appendix A. 4. Secondly, the convergence of the perturbative expansion at the low scale 

~ A may be improved by borrowing the idea of "contour resummation" developed in [68]. 
Since the shape-function integrals can be written as contour integrals in the complex plane 
along a circle with radius A, in may be more appropriate to use a contour-weighted coupling 
constant rather than the naive coupling q;s(A). Exploring the numerical impact of these two 
proposals is left for future work. 
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Appendices 

A.l Anomalous dimensions and RG functions 

The exact solutions (14) to the RG equations in (12) can be evaluated perturbatively by 
expanding the anomalous dimensions and j3 function, 



rcusp('^s 



I3{as) 



-2a, 



471/ 
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/50 7^+/5l 7^ 



47r 



47r 



(59) 



and similarly for the remaining anomalous dimensions. We work consistently at next-to- 
leading order in RG-improvcd perturbation theory, keeping terms through order ag in the 
final expressions for the Sudakov exponent S and the functions ar, a-y, and ay. For ar one 
obtains the standard expression 
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The result for the Sudakov factor S is more complicated, as it is necessary to include terms 
of next-to-next-to-leading logarithmic order. We obtain 
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where r = as{ii)/as{i'). Whereas the two-loop anomalous dimensions and /3 function are 
required in (60), the expression for S also involves the three-loop coefficients r2 and /32- 

The perturbative expansion of the QCD /? function to three-loop order is [69] (all results 
refer to the MS renormalization scheme) 
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where the numerical values refer to Ai'c = 3 and n/ = 4. The two-loop coefficient of the cusp 
anomalous dimension has been known for a long time [33]. However, its three-loop coefficient 
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Figure 3: Sudakov exponents S{mb, /i) (black) and S{1 GeV, /i) (gray) at next- 
to-next-to- leading order (solid), next-to-leading order (dashed), and leading 
order (dash-dotted) . The sohd and dashed curves are nearly indistinguishable. 

has only been calculated very recently by Moch et al. [34]. This is a lucky coincidence, because 
that calculation was done in a context not related to heavy-quark physics. The results are 

ro = 
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Although the two- and three- loop coefficients of the /3 function and cusp anomalous anomalous 
dimension are large, the perturbative expansion of the Sudakov exponent is extremely well 
behaved. This is illustrated in Figure 3, which shows the Sudakov exponents S{iJ,h, fi) and 
S{i-io, ij) for Hh = nib = 4.7 GeV and Ho = I GeV as a function of /i. 

The two-loop coefficient of the anomalous dimension 7 entering the shape-function evolu- 
tion kernel in (17) has been calculated in [64]. We have found some mistakes in the translation 
of the results for the two-loop graphs into the expression for the anomalous dimension. The 
corrected result is [70] 

7o = -2C^--|, 



71 = -8C, 



37 TT^ 9 , 

h - Cs 

108 144 4 



-66.7531 + 4k , (64) 



where k, = under the assumption that the two-loop diagrams themselves were evaluated 
correctly in [64] . However, there is reason to believe that there might be an additional error in 
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this paper, giving rise to a non-zero value k — 4/3 [65] (see also Appendix A. 2 below), which 
we adopt in our numerical analysis. 

The two-loop anomalous dimension for the leading-order SCET current operator in (13) 
has not yet been computed directly. An analysis is in progress and has already led to a 
prediction for the terms of order Cpfif [37]. The remaining terms can be deduced by noting 
that the difference 7"^ = 7' — 7 is the non-cusp part of the anomalous dimension of the jet 
function [12], which is related to the familiar jet function from deep-inelastic scattering. We 
find 



(-^ + 27r2-24C3)c^+(- 



3 155 227r2 
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(65) 



where the terms in the first line in the expression for 7j^ are taken from [35] , while the remainder 
in the second line is due to the non-trivial normalization of the SCET jet function in (9). 
Combining the results (64) and (65), we obtain 
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Only the term proportional to n/ in 7^ has so far been checked by a direct calculation in SCET 
[37]. 

A. 2 Pert ur bat ive expansion 

In this work, we have presented for the first time the complete RG-improved expression for the 
B Xs7 decay rate, in which all logarithms In 5 (with 5 ~ A/m^) are resummed at next-to- 
next-to-leading logarithmic order. In order to simplify the comparison of our result with those 
of other authors, we expand it in fixed-order perturbation theory and list the resulting terms 
at order ag and a^. It suffices to focus on the perturbative correction to the term multiplying 
the product m^m^(///i) |C7-y(//fe)p, where is the pole mass. We find 
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where 
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These expressions are independent of the matching scales /ij and /iq, and they have the correct 
dependence on The constant ko can only be obtained from a complete calculation of 0{al) 
corrections to the decay rate. 

Resummed expressions for the B — > Xg^ photon spectrum with next-to-leading logarithmic 
accuracy have been reported in [71] and [36]. In the first paper, expressions for the coefficients 
^4 and ^3 are derived, which agree with our findings. In [36], Gardi has obtained a result for 
the photon spectrum from which all four coefficients ki can be extracted. By matching his 
result for ki with ours, we conclude that k — 4/3. 



A. 3 Kinematic power corrections for the average photon energy 

The functions dij{S) entering the expression for the average photon energy in (49) are given 

by 
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where xs — max(a;, 1 — 5), z — [mc/mi^Y, and the function G{t) has been defined in (38). 
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A. 4 Comment on BLM scale setting 



Here we illustrate the simple fact that the BLM scale-setting procedure [66] for multi-scale 
problems is more complicated than in the familiar case with a single scale. Let us, for simplicity, 
ignore RG resummation effects due to anomalous dimensions and consider a physical quantity 
A, whose perturbative expansion is given by the product of two perturbative series at scales 
M and m, with M > m. We write 

1 + ci a{M) + (2/3oC2 + 4) a^{M) + ...]\l + di a{m) + (2/30^2 + 4) + 

(70) 

where a = as/{4:'K), and cf\ df' are numerical coefficients. The BLM scales of the two series 
are 

/^b!m = M e-^/^^ , /iL°LM = m e''^"^ . (71) 

They are determined so as to absorb the 0{a^g) terms multiplying /5o into the running coupling 
constants. Adopting the BLM philosophy, we would conclude that perturbation theory is well 
behaved as long as both /^blm ^^<i A*blm in the perturbative regime. 

Imagine now that we compute A in fixed-order perturbation theory using a single coupling 
constant as{^J)■ We would obtain 



^ = 1 + (ci + di) a{ij) + 
and the associated BLM scale would be 



r^n ( 7 1-^ , ^ ^\ II II 1 \ 
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a2(/x) + ... , 
(72) 



^BiM = M[-) exp(^--^j. (73) 

Obviously, equation (72) does not provide the same information as (70), and in particular it 
does not allow us to compute the BLM scales in (71). To this end we would need C2 and ^2 

separately, not just their sum. 

It is instructive to look at a couple of examples, where the conclusions derived from (72) 
would differ strongly from those derived from (70). Consider, for instance, a situation where 
(ci + di) is accidentally small. Then the BLM scale (73) is either very small or very large, 
whereas the BLM scales in (71) could be close to the scales M and m, respectively. BLM scale 
setting based on the fixed-order calculation would then be totally misleading. Next, consider 
the case where the coefficient of (3o in (72) is small, for instance by a particular choice of /x. 
The fixed-order calculation would lead us to conclude that BLM-type corrections are small, 
whereas the BLM-type terms in (70) could still be large. Finally, consider the (somewhat 
pathological) example where di ~ ci and ^2 — — C2 with 7 = ^2/(^1 > 0. Then the "physical" 
BLM scales are /Ublm — ^'^M and /UgLM — 7 large, perturbation theory may be in 

trouble, since /isLM may no longer be in the perturbative regime. Nevertheless, the "average" 
BLM scale /^blm — V Mm is large, and the fixed-order calculation would thus indicate a 
well-behaved perturbative expansion. 
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